Three-body recombination in a three-state Fermi gas with widely tunable interactions 
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We investigate the stability of a three spin state mixture of ultracold fermionic ^Li atoms over 
a range of magnetic fields encompassing three Feshbach resonances. For most field values, we 
attribute decay of the atomic population to three-body processes involving one atom from each 
spin state and find that the three-body loss coefficient varies by over four orders of magnitude. We 
observe high stability when at least two of the three scattering lengths are small, rapid loss near the 
Feshbach resonances, and two unexpected resonant loss features. At our highest fields, where all 
pairwise scattering lengths are approaching at — — 2140ao, we measure a three-body loss coefficient 
La ~ 5 X 10~^^ cm^/s and a trend toward lower decay rates for higher fields indicating that future 
studies of color superfluidity and trion formation in a SU(3) symmetric Fermi gas may be feasible. 

PACS numbers: 67.85.Lm, 34.50.-s, 67.85.-d, 03.75. Ss 



Multi-component Fermi gases with tunable interac- 
tions are exceptionally well suited to the study of 
few- and many-body quantum physics. Ultracold two- 
state Fermi gases near a Feshbach resonance have been 
used to characterize the crossover from Bardeen-Cooper- 
SchreifFer (BCS) superfluidity to Bose-Einstein conden- 
sation (BEC) of diatomic molecules [J. Further, a nor- 
mal to superfluid transition and phase separation have 
been observed in imbalanced two-state spin mixtures . 
The stability of the two-state mixtures against two- and 
three-body loss processes was critically important to the 
success of these experiments. 

Recently, there has been considerable interest in the 
study of three-state Fermi gases with tunable interac- 
tions |3j S S 0j Si ■ Whereas three-body interactions 
in ultracold two-state Fermi gases are suppressed by the 
exclusion principle, the addition of a third spin compo- 
nent allows for the study of three-body phenomena such 
as the Efimov effect Q ■ Further, with this additional spin 
state there may be competition between multiple pairing 
states and trion formation 0, Hi- For example, if pair- 
wise interactions are all attractive and of equal magni- 
tude, the system is expected to exhibit a novel superfluid 
phase analogous to color superconductivity in quantum 
chromodynamics (QCD) A quantum phase transition 
to a Fermi liquid of trimers (analogous to baryons) may 
be observed if the ratio of the interaction energy to the 
kinetic energy is increased (e.g. by an optical lattice) [^. 

Future studies of the above phenomena critically de- 
pend on the magnitude of two- and three-body loss rates, 
particularly when two or more scattering lengths are reso- 
nantly enhanced. Two- and three-body loss and heating 
processes can impose stringent limits on the maximum 
achievable phase space density. Further, the unambigu- 
ous observation of three-body resonances requires negli- 
gible two-body loss 

Most investigations of Fermi gases with tunable inter- 
actions have focussed on two-state mixtures. A small 
admixture of a third spin component has been used 



for thermometry [lO[ and recently the rapid decay of 
a three-state Fermi gas near a Feshbach resonance was 
noted during an investigation of Cooper pair size by 
radio-frequency (RF) spectroscopy in a two-state Fermi 
gas [m . Very recently our group reported measurements 
of three-body loss in thermal and degenerate three-state 
Fermi gases of ^Li for magnetic fields between 400 and 
960 G which included three Feshbach resonances [l^ . 
Independently, Ottenstein et al. observed loss in a three- 
state ^Li gas for field values between - 750 G (includ- 
ing a single Feshbach resonance) and reported three-body 
loss rate coeflacients for field values up to 600 G (exclud- 
ing any pairwise scattering resonances) [l3| . 

In this Letter, we report the first measurements of 
three-body loss coefficients for a three-state Fermi gas 
over a range of magnetic fields where pairwise interac- 
tions are resonantly enhanced. This enhancement is due 
to three overlapping Feshbach resonances and a zero- 
energy resonance in the ®Li triplet molecular potential. 
Further, we demonstrate that all two-state mixtures of 
the three lowest energy hyperfine states of ®Li are sta- 
ble against two-body loss processes for almost the entire 
range of field values between 15 and 953 G, making this 
three-state mixture well suited to studies of three-body 
physics. We observe a narrow loss feature at 127 G [1J| 
which may be due to a three-body resonance and an addi- 
tional state-dependent loss feature at 504 G. Finally, we 
measure the three-body loss rate coefficient at high fields 
where all three scattering lengths are asymptoting to the 
triplet scattering length at = — 2140ao. This measure- 
ment and an observed trend toward lower loss at higher 
fields has important implications for realizing cold atom 
analogs of color superfluid and baryon phases in QCD. 

We study a ^Li Fermi gas with equal populations in the 
three lowest energy hyperfine states. At zero field, the 
three states correspond to |1) = 1^,-1-^), |2) = |^,— 5) 
and |3) = in the |/, m/) basis. For fields above 

~ 200 G, these states become increasingly electron-spin 
polarized and primarily differ by their nuclear spin pro- 
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FIG. 1: (color online) The s-wave scattering lengths as a func- 
tion of magnetic field for interactions between the three lowest 
energy hyperfine ground states of ^Li |2) and |3)) [l5l |. 
Feshbach resonances occur at 690, 811, and 834 G. 

jection. The pairwisc s-wave scattering lengths between 
these states (ai2, 023, and 013) exhibit three broad over- 
lapping Feshbach resonances as shown in Fig.[T][l^. For 
very large magnetic fields, each of the pairwise scattering 
lengths asymptote to the triplet scattering length for ^Li, 
at = — 2140ao. Therefore, in the limit of large magnetic 
fields, the interactions are SU(3) symmetric [5|. 

Two-body spin-fiip processes are expected to be espe- 
cially small for this ultracold ^Li mixture. Spin-exchange 
collisions are energetically forbidden in a magnetic field. 
Dipolar relaxation by electron spin-flip is suppressed at 
high field as the gas becomes increasingly electron-spin 
polarized in the lowest energy electron spin state. 

We can produce a degenerate Fermi gas (DFG) with 
an equal mixture of atoms in states |1) and |2) once ev- 
ery 5 seconds by evaporatively cooling this mixture in 
a red-detuned optical dipole trap [l^. During the first 
second, ~ 10* ^Li atoms from a Zeeman-slowed atomic 
beam are collected in a magneto-optical trap (MOT). 
A crossed optical dipole trap which overlaps the MOT 
is then turned on and optimally loaded by tuning the 
MOT laser beams ~ 6 MHz below resonance and reduc- 
ing their intensity for 7 ms. The atoms are then optically 
pumped into states |1) and |2) prior to extinguishing the 
MOT laser beams and field gradient. The two beams 
which form the optical trap derive from a linearly polar- 
ized multi- longitudinal- mode 110 W fiber laser operating 
at 1064 nm. The beams nearly co-propagate in the verti- 
cal direction (y), have orthogonal polarizations and cross 
at an angle ~ 11°. The beams are elliptical with calcu- 
lated waist radii ^ 30/im and ^ 100/im at the point 
of intersection. The maximum trap depth per beam is 
~ 1 niK allowing ~ 5 x 10^ atoms to be initially loaded. 
We apply a bias field and a noisy RF pulse to create a 
50-50 mixture of atoms in states |1) and |2) [lit . 

Forced evaporation of the atoms occurs at a field 
of 330 G where the two-state scattering length ai2 — 
— 280ao. We lower the depth of the optical trap Uq by a 
factor of 107 over 3.6 seconds to obtain our final temper- 
ature. For the work presented here we study a gas at a 
temperature T > 0.5 Tp {Tp is the Fermi temperature) 
so that it is appropriate to treat the cloud as a thermal 
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FIG. 2: (color online) The population in states |2) and |3) 
measured by absorption imaging in the Paschen-Back regime 
at a bias field of 568 G [TH]. A a~ polarized probe beam drives 
nij = — I — * "T-j' = — I transitions which, due to the hy- 
perfine interaction in the ground state, are spectroscopically 
resolvable. (Top right) The on-resonance absorption images 
show the density profile for the three states following a 900 /is 
time-of-flight expansion. (Bottom right) The total absorption 
vs. probe detuning (relative to the [/ = ^,mf = — ^) 
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) transition frequency in zero field). 



gas in the analysis. To suppress further loss by evapora- 
tion for the remainder of the experiment, Uq is increased 
by a factor of 4 increasing the ratio of Uo/kBT by a fac- 
tor ~ 2. The final oscillation frequencies of the trap are 
ly^ = 3.84 kHz, i^y = 106 Hz and = 965 Hz [l3| with a 
final trap depth per beam ~ 40 /xK. The total number 
of atoms ~ 3.6 x 10'' in a balanced two-state mixture 
at T ~ 1.9 /xK for which Tp ~ 3.7/iK. The background 
limited 1/e-lifetime ~ 30 s. 

To create an incoherent three-state mixture we first 
increase the strength of the magnetic field in 10 ms to 
568 G where the mixture is stable. We then apply a 
noisy RF magnetic field with two frequencies centered on 
the |1) — 1 2) and |2) — |3) transitions (each with a power 
corresponding to an on-resonance Rabi frequency ^ 
27r X 18kHz and both broadened to a width of 1 MHz). 
A magnetic field gradient of ~ 1 G/cm in the direction of 
the bias field (z) is simultaneously applied to destroy any 
internal state coherence in the sample. The noisy RF field 
remains on for 50 ms to create an equal mixture of iV ~ 
1.2 X 10^ atoms in each state. At this point, T = 1.9 
with Tp = 3.2 /xK and the peak density in a single spin 
state is 5.5x 10^^ atoms/cm'^. We measure the population 
and temperature of atoms in each state separately by 
spectroscopically-resolved absorption imaging (Fig. [2]). 

To investigate the magnetic field dependence of decay 
in this three state mixture, we compared the number of 
atoms remaining in each spin state after spending two 
different times, 1 ms and 201 ms, at a particular field 
of interest B (see Fig. El^a)). For each experimental cy- 
cle the magnetic field was ramped in 10 ms from 568 G 
(where the mixture was first created) to the field of inter- 
est where it was held constant for either 1 ms or 201 ms. 
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FIG. 3: (a) Number remaining in each of the three spin states 
after spending 201 ms at the field of interest B (normalized 
to the number remaining after spending 1 ms at B). (b) 
Magnetic field dependence of decay in individual two-state 
mixtures. ♦ (•) indicates the normalized number remaining 
in state |2) (|3)) of a |2> - |3) (|1) - |3>) mixture after 201 ms. 
The thick lines guide the eye. The representative error bars in- 
dicate the standard deviation in the mean. Absorption images 
are taken at 953 G to ensure that any weakly-bound molecules 
which remain trapped are dissociated and measured. 

The field was then ramped in 10 ms to 953 G and held for 
20 ms before absorption imaging at 953 G to measure the 
atom number (A^ims or iV2oims)- The 10 ms field ramp to 
953 G (where 012,023,013 < 0) ensured that any atoms 
which had formed weakly bound molecules in the vibra- 
tional state associated with the Feshbach resonances (but 
remained trapped) would be dissociated and measured. 
For each field value, A^20ims is normalized by -/Vims to cor- 
rect for the fraction of atoms lost during the field sweeps, 
both to and from the field of interest B, which exhibited 
a slowly varying dependence on B. See Rcf. [2l| for the 
data before normalization. 

To confirm that two-body loss does not contribute sig- 
nificantly to the measurements described above, we mea- 
sured the ratio of atoms remaining after evolution times 
of 201 ms and 1 ms at fixed magnetic field values for 
each of the three possible two-state mixtures. Figure[3l^b) 
shows the ratio A*2oims/-^ims for atoms in state |3) (state 
|2)) of a |1)-|3) (|2) -|3)) mixture. Data for the |1)-|2) 
mixture (not shown) is very similar to that of the |2) — |3) 
mixture (e.g. Ref. [19[). With the exception of loss fea- 
tures between 600 and 750 G, each two-state mixture was 
stable. The loss features are due to three-body recombi- 
nation to bound molecular states associated with the Fes- 
hbach resonance. Since the two-state mixtures are stable 
for fields < 600 G and > 750 G, loss features observed in 
three-state mixtures at these fields are due to three-body 



events involving one atom from each spin state. 

In Fig. El^a), the broad dominant loss feature centered 
at 720 G occurs in the vicinity of three overlapping in- 
terspecies Feshbach resonances. Significant loss due to 
three-body recombination is expected near these reso- 
nances since recombination events involving one atom in 
each spin state are not suppressed by the exclusion prin- 
ciple and one expects a significant increase in the event 
rate when two or more scattering lengths are resonantly 
enhanced [2^. Similarly, high stability at zero field and 
near the zero crossings of the Feshbach resonances is not 
surprising since at least two of the scattering lengths are 
small at these fields. At high field, the stability increases 
relative to that at 720 G. 

Unexpected resonant loss features are observed at 127 
and 504 G, where the two-body scattering lengths are 
not predicted to exhibit any resonances. We have ob- 
served these features at hi gher temperatures and in a 
different trap configuration [21j. A possible explanation 
for the feature at 127 G is that the binding energy of 
a three-body bound state crosses through zero near this 
field value. At 504 G, we observe differing degrees of 
enhanced loss for the three states leading to population 
imbalance. This resonance was not observed in Ref. [Tst - 

The feature at 228 G in Fig. [3] is consistent with our 
estimated location of a |1) — |3) p-wave Feshbach reso- 
nance which would yield the observed preferential loss of 
atoms from states |1) and |3) [i^l- We also observed a 
very narrow but inconsistent loss feature at 259 G (not 
shown). The modest increase in stability observed near 
200 G did not appear in our other experimental data [2l[ . 

A second set of experiments measured three-body loss 
coefficients at fixed magnetic fields where the decaying 
populations remained balanced. In these experiments, 
the magnetic field is swept to the field of interest B in 
10 ms after creation of the three-state mixture and the 
number N{t) and temperature T{t) of atoms in state |3) 
are then measured at the field B for various delay times 
t. To extract the three-body loss coefficients, we fit this 
data assuming one- and three-body loss and heating due 
to three-body recombination. 

When three-body recombination involves one atom 
from each spin state, the number of trapped atoms in 
each of the equally populated spin states N{t) evolves 
according to iV = ~LiN — L^{'n?)N. Here, {n^) is the 
average value of the squared density per spin state and 
L3 (Li) is the three-body (one-body) atom- loss rate co- 
efficient. The density distribution in the trap is well de- 
scribed by a thermal distribution so that 
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where 7 = L3(m(I'^/27rfcB)'^/v27 and uj ~ [ujxUJyUOzY 
Following Ref. [2^, we model the temperature increase 
in the gas as arising from "anti-evaporation" and recom- 
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FIG. 4: Three-body loss rate coefficient, is, at various fields, 
bination heating. The temperature evolves as 
dT {T + n) 

where fesT^ is the average energy deposited in the gas 
per recombination event. 

For a given parameter set {N{0), T{0), Li, and 
T/j), Eqns.[T]and [2]can be numerically integrated to give 
N{t) and T{t). The best-fit parameters for a data set at 
a given field are determined by minimizing x^. The best- 
fit L3 values at various fields are displayed in Fig. 21 For 
these fits, Tfi varied from 1.5 — 4 fiK and Li varied from 
0.1 - 0.33 s^i. The error bars in Fig. H indicate the 68% 
confidence interval due to statistical uncertainty in all fit 
parameters. The error bars do not include a systematic 
uncertainty of ±70% arising from our uncertainty in atom 
number (±30%) and trap frequencies (±5%). 

As shown in Fig. [H L3 varies by over four orders of 
magnitude. In the range 600 G < B < 750 G, the L3 val- 
ues (o) may be overestimated since the measured decay 
includes loss events we observe in two-state mixtures at 
these field values. Near and 568 G, three-body recom- 
bination rates are relatively small and stable three-state 
Fermi gases can be created. For L3 = 10~^^cm^/s, a 
gas with a density n — 10^^ cm~'^ per spin state has 
a 1/e-lifetime 10 seconds. Consistent with the data 
shown in Fig. El^a), a resonant peak in L3 is observed 
at i? ~ 127 G. Wc do not report a value for L3 at the 
location of the 504 G resonant feature (Fig. [31(a)) since 
a population imbalance develops there. The highest L3 
value we report (lO^^^cm^/s) occurs at 775 G and is 
likely unitarity limited (for fields > 600 G, T ~ 6 /xK 
due to heating). We observe that L3 decreases by a fac- 
tor of ~ 2.5 as the field is increased from 775 to 953 G, 
the highest field we currently access. This trend suggests 
that L3 may decrease further in the high field limit where 
ai2, 023, ai3 ~* -2140ao. 

In conclusion, we can create long-lived three-state 
DFGs of ^Li atoms for fields near 568 G. The minimal 
two-body loss we observe makes this system promising 
for future studies of Efimov trimers. Indeed, the loss 
resonances at 127 and 504 G may be occurring at field 



values where trimer states cross the dissociation thresh- 
old. The recombination rates at high fields need not pre- 
clude experiments designed to study color superfluidity 
and trion formation For instance, using low density 
gases (n ~ 5 X 10^° cm"-^) and large period {d ~ 2 ^m) 
optical lattices, long lifetimes (> 0.1s) and strong inter- 
actions can simultaneously be achieved. 
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